Atypical BCS-BEC crossover induced by quantum-size effects 
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Quantum-size oscillations of the basic physical characteristics of a confined fermionic conden- 
sate are a well-known phenomenon. Its conventional understanding is based on the single-particle 
physics, whereby the oscillations follow variations in the single-particle density of states driven by 
the size quantization. Here we present a study of a cigar-shaped ultracold superfluid Fermi gas, 
which demonstrates an important many-body aspect of the quantum-size coherent effects, over- 
looked previously. The many-body physics is revealed here in the atypical crossover from the 
Bardeen-Cooper-Schrieffer (BCS) superfluid to the Bose-Einstein condensate (BEC) induced by the 
size quantization of the particle motion. The single-particle energy spectrum for the transverse di- 
mensions is tightly bound, whereas for the longitudinal direction it resembles a quasi-free dispersion. 
This results in the formation of a series of single-particle subbands (shells) so that the aggregate 
fermionic condensate becomes a coherent mixture of subband condensates. Each time when the 
lower edge of a subband crosses the chemical potential, the BCS-BEC crossover is approached in 
this subband, and the aggregate condensate contains both BCS and BEC-like components. 

PACS numbers: 67.85.Lm 



I. INTRODUCTION 



A crossover from the wcakly-intcracting BCS super- 
fluid of Cooper pairs to a BEC of tightly-bound molecule- 
like pairs m y| is one of the most important phenomena 
in the physics of fermionic condensates (see also the re- 
view [3|). Although the BCS-BEC crossover vifas origi- 
nally discussed in the context of semiconductor materials 
in the presence of superconducting correlations [l[ (for 
recent activity on this subject, see |j, Q), it was first 
demonstrated in experiments with ultracold superfluid 
Fermi gases Q. In the standard scenario the BCS-BEC 
crossover for ultracold fermions is achieved via a Fesh- 
bach resonance in the particle scattering. Here we report 
a different mechanism to achieve the BCS-BEC crossover, 
inducing it by quantum-size (QS) effects in a cigar-shaped 
superfluid Fermi gas. 

The QS effects in superfluid or superconducting sys- 
tems have their origin in the geometric quantization of 
the single-particle motion and reveal themselves in the 
oscillations of the basic condensate properties. Such os- 
cillations have been theoretically investigated in many 
different systems e.g., ultrathin films [6|; quantum striped 
superconductors and superconducting heterostructures 
at the atomic limit [3, Q; superconducting metallic 
nanowires [9|; and a pancake-shaped superfluid Fermi 
gas [lO| . Recently, atomically uniform Pb nanofllms were 
fabricated, which resulted in the flrst experimental ob- 
servation of the QS oscillations in the critical temper- 
ature |ll| . An interesting experimental study of ^Li 
Fermi gas in a pancake-shaped trap was also recently 
reported [I4I where the effects of the transverse quantiza- 
tion on the aspect ratio of the atomic cloud were demon- 
strated. This opens new prospects for the study of the 



QS oscillations of the properties of fermionic condensates 
in ultracold Fermi gases with tunable confinement pa- 
rameters. Conventional understanding of the QS oscilla- 
tions follows from the fact that the single-particle energy 
spectrum for the quantum-confined dimensions is tightly 
bound, whereas along the other dimensions it resembles 
a quasi-free dispersion. This results in the formation of a 
series of single-particle subbands. The lower edge (bot- 
tom) of such a subband coincides with the corresponding 
discrete single-particle level associated with the confined 
dimensions. When the energy spacing between subbands 
is systematically decreased, e.g., by decreasing the rele- 
vant trapping frequency for superfluid fermionic atoms or 
the thickness of an atomically-uniform metallic nanofilm, 
the subband lower edges sequentially cross the chemical 
potential /x. Each time this happens, the single-particle 
density of states (DOS) at fi increases, leading to a higher 
critical temperature, larger excitation gap, etc., which is 
referred to as the shape or QS resonance [a, H • 

In the present work, based on a study of a cigar-shaped 
ultracold superfluid Fermi gas, we demonstrate that the 
QS coherent effects cannot be fully understood in terms 
of the single-particle physics. The total or aggregate 
fermionic condensate in the system of interest is a coher- 
ent mixture of the subband components (condensates), 
and each component undergoes a BCS-BEC crossover 
when the lower edge of the corresponding single-particle 
subband crosses fi. As a result of such an atypical BCS- 
BEC crossover, the total condensate is a coherent mix- 
ture of both the BCS and BEC-like components, which 
is most pronounced at the shape resonances. 

This paper is organized as follows. In Sec.|IT]we outline 
the relevant formalism for calculating the "wave func- 
tion" of a condensed fermionic pair in the cigar-shaped 
superfluid Fermi gas. In Sec. IIIII we analyze numeri- 



cal results for the fermionic-pair "wave function" with 
the focus on the atypical BCS-BEC crossover induced by 
the QS effects as dependent on the perpendicular (trans- 
verse) trapping frequency. In Sec. llVl we discuss techni- 
cal details of the approximations used in our approach 
and show that the results here reported are not sensitive 
to these approximations. Our conclusions are given in 
Sec.lVl 



II. FORMALISM 

Our analysis is based on a numerical solution of the 
Bogoliubov-de Gcnnes (BdG) equations for ^Li atoms 
trapped by a harmonic axially symmetric potential 
[/(r) = M(w5^yO^ -|-a;| 2:^)/2, with r = {p, (j), z} cylindrical 
coordinates and uj\\ ^ uj±. Calculations are performed 
for zero temperature as the BdG equations are appro- 
priate to describe the BCS-BEC crossover in spatially 
nonuniform fermionic systems at nearly zero tempera- 
tures. In particular, it has been shown in [13| that the 
BdG equations reproduce the Gross-Pitaevskii equation 
for the condensate wave function on the BEG side of the 
crossover. 

To clearly display the underlying physics, we employ 
the Anderson semi-analytical approximation [IJ] accord- 
ing to which the spatial dependence of the particle-like 
w^(r) and hole- like w,y(r) wave functions appearing in the 
BdG equations is chosen to be proportional to the cor- 
responding single-particle wave function '^jy(r) (see the 
discussion in Sec.lIVp. i.e.. 



My(r) = U^ip^ir), Wy(r) = V^ip^{r), 



(1) 



where ly — {j, n, m} are the three quantum numbers as- 
sociated with the longitudinal motion along the z axis, 
the radial motion and the angular momentum, respec- 
tively. iy9,y(r) = •dnm{p,4>)Xj{^) IS the product of the 
eigenfunctions of the 2D and ID harmonic oscillators. 
Inserting Eq. ([T]) into the BdG equations, one gets a sys- 
tem of two linear equations for the coefficients Ui, and 
Vu (chosen real). A nontrivial solution exists when the 
corresponding determinant is equal to zero, which gives 
the quasiparticle energy E^, = a/A^ + AJ, where \i, — 
hjjj±{l + 2n-\- \m\)-\-fvjj\\{j + \/2) — ^ is the single-particle 
energy measured from the chemical potential, and IS.^ is 
the corresponding pairing energy. Then, together with 
the normalization condition Z-/^ -f V^ ~ 1 , the BdG equa- 
tions yield Ul ^ {1 + K/E,)/2, Vl = (l - A,/S,)/2. 
These expressions for lAu and Vi, make it possible to 
find the BCS-like self-consistency equation given by (at 
T = 0) 



= ^E^-A.,(^-^), 



(2) 



with the interaction matrix Vyi^i = 

(7/d'^r|(/5j^(r)|-^|(/5i^/(r)p and g the coupling constant. 
The second term in brackets of Eq. ^ eliminates the 



ultraviolet divergence: this is a convenient simplification 
that, to first approximation, models the rigorous reg- 
ularization for a spatially nonuniform system reported 
in [131 (see the discussion in Sec. lIVp . We avoid the 
ultra-confinement regime where the effective dimension- 
ality of the system reduces and the particle scattering 
becomes different from the 3D case (see, e.g., Q]). Here 
we take /i > 2fiLo±^ and, in addition, the absolute value 
of the s-wavc scattering length a {a < 0) is chosen 
smaller than ^|[,^j_, where la = \/h/{MuJa)- With this 
choice the interatomic collisions can be regarded as a 3D 
process (see the discussion in Sec. IIVP for which we use 
the standard expression of the pseudopotential theory 
g = 47rfi,^|a|/M, with M the atomic mass. 

Equation ([2|) can be viewed as a system consisting of 
multiple condensates with the pairing gaps Ai/ coupled 
through the interaction matrix Vi/^'. Furthermore, for 
the cigar-shaped trap the interlevel energy spacing cor- 
responding to the quantization in the z direction is suf- 
ficiently small so that the single-particle spectrum can 
be viewed as a sequence of the subbands (n, m). Indeed, 
differences between Aj^'s within the same subband are 
almost insignificant (and disappear in the limit lu —^ cxj). 
It is therefore useful to distinguish separate subband 
contributions to the system characteristics, i.e., to treat 
the system as a coherent mixture of multiple subband- 
dependent pair condensates. 

By solving Eq. ^ we obtain the set of A^. To probe 
the spatial pairing correlations (the main point of our 
study), the anomalous correlation function ^(r,r') = 
('0^(r)f/;|(r')) need to be calculated. Following the orig- 
inal works of Gor'kov[16| and Bogoliubov l3l, it can be 
viewed as the wave function of a condensed fermionic 
pair. Using our subband-based classification, we can rep- 
resent \I'(r,r') (using the Bogoliubov canonical transfor- 
mation) as a sum over the relevant subbands 



vl/(r,r') = 5] *„™(r,r'). 



(3) 



where *„m(r,r') = dnm{p,4>Wnm{p' ^4>') ipnm{z,z') and 
at zero temperature 

'fpn7n{z, Z') =- ^ Xj{z) X]{z') ^nmj 



1 



1 



^rnnj ■^iimj 



(4) 



Due to pairing, \I'(r, r') is localized as a function of the 
longitudinal relative coordinate z — z' (in the x, y plane it 
is confined by the trapping potential) and the character- 
istic localization length, i.e., the longitudinal fermionic- 
pair size, is calculated as 



^0= (aA-1 /"d3rd3r'|*(r,r')P(z-z')') 



1/2 



(5) 



where M = J'd'^rd^r'|^(r,r')p is the normalization fac- 
tor. Similarly, one can define a subband-dependcnt 



fermionic-pair size ^^"'"^ given by Eq. ([5]), with ^(r, r') 
and A/" replaced by \l/„m(r, r') and the corresponding nor- 
mahzation factor Mnm ■ We note that while other defini- 
tions of the condensed-pair size are possible, the resulting 
expressions differ only by a constant factor [1^, [I^l ■ 



III. ATYPICAL BCS-BEC CROSSOVER 

As an illustration we consider a mixture of ^Li 
fermionic atoms with two interacting spin states, 
\F,mF) = |l/2,l/2) and |f/2,-l/2). The correspond- 
ing scattering length a can be significantly modified 
through a broad Feshbach resonance with experimen- 
tally reported values on the BCS side from — 250nm to 
— 100 nm [201 ■ For our numerical calculations we take a = 
-140 nm and -180 nm. We choose a;||/(27r) = 240 Hz, 
/i = 100fta;||, and uj±_ is assumed to be variable leading to 
tunable QS effects. (We found it convenient to present 
our results versus the ratio s = ^/hu)±). For the above 
choice of the parameters the particle density in the cen- 
ter of the trap is found to be about 10^^-10^"^ cm~^, with 
kp\a\ = 0.7-0.9. This is consistent with most experi- 
ments where the particle density is reported in the range 
10i2_ioi5cm-3 and kpH ~ 1 [3. The characteristics 
of our trapping potential are chosen similar to those re- 
alized in recent experiments with quasi-lD Fermi gases. 
In particular, one can compare a;n^'^/(27r) = 200 Hz and 

the ratio cjjj^'^/wl]^'^ = 0.001 reported in [2l| with our 
values 240 Hz and 0.02-0.04. Note that u}\\/u}± is larger 
in our calculations because, as already mentioned above, 
we avoid the ultra-confinement regime for the transverse 
motion of atoms that is realized in J2l| and where only 
one transverse level (i.e., one subband) is occupied. 



-140 nm 



A. Oscillations of longitudinal fermionic-pair size 

To reveal the many-body aspect of the QS effects, we 
first need to consider how A^'s change with s = ^/hu)±. 
Figures [TJa) and (c) show (A) i.e., the pairing gap aver- 
aged over the Fermi-surface, as a function of s calculated 
for a = — 140nm and — 180nm, respectively. As can be 
seen, (A) increases in the vicinity of integer values of s. 
This condition for the developing of a shape resonance [6[ 
is satisfied when the bottom of a subband, referred to as 
the resonant subband, approaches the chemical potential 
and the DOS increases. As a result, (A) exhibits QS os- 
cillations with changing uj± similar to those reported for 
a pancake-shaped superfluid Fermi gas |1G |. 

Now, based on the data for (A), we demonstrate that 
the conventional single-particle picture of the QS oscil- 
lations fails to explain the corresponding changes in the 
two-particle characteristics. This is illustrated by the re- 
sults for the longitudinal fermionic-pair size ^g shown in 
Figs. [IJb) and (d). Similar to the averaged pairing gap, 
^0 exhibits remarkable QS oscillations. However, using 
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FIG. 1: (Color online) Panels (a) and (c) demonstrate the 
QS oscillations of the averaged pairing gap (A) for scatter- 
ing lengths a = — 140nin and — 180nm, respectively. Pan- 
els (b) and (d) show the corresponding changes in the longi- 
tudinal fermionic-pair size ^o [see Eq. ([5}] and the subband 
fermionic-pair size q^"^> for subbands (n, m) = (0, 0), (0, ±1) 
and (0,±2). 



the standard BCS estimate ^o oc hvp/{A)^ with vp the 
Fermi velocity at the center of the trap, and taking into 
account the variations of (A) in Fig. [U we obtain for ^q 
a decrease by a factor of 1.6 when s increases from 1.8 to 
2.1 at a = —140 nm. Note that the corresponding change 
in vp is negligible, see the data for single-particle density 
discussed in Sec. IIII CI This is a considerable underesti- 
mation of the numerical results given by Fig.[ljb), where 
Co|s=i.8/'?o|s=2.i ~ 3. A similar discrepancy is found for 
a = —180 nm. 

A detailed analysis shows that this discrepancy is re- 
lated to a significant redistribution of the fermionic con- 
densate over the available subbands. This is seen from 
Fig. nib), where ^o is compared with ^^"™^ for subbands 
(n, m) = (0, 0), (0, ±1), (0, ±2). When s < 1.9, the main 
contribution to the pair condensate comes from subband 
{n,m) = (0,0), i.e., \1/ « ^o,0: and we obtain Co ~ Co ■ 
As the system goes through the resonance that develops 

at s = 2, Co drops and approaches Co -In this case, 
two resonant subbands (n, m) = (0, ±1) make the largest 
contribution to the total condensate i.e., about 70% at 
s = 2.1. At the next resonance, s = 3, Co decreases again 
and approaches Co : which points to the enhancement 
of the contribution ^o,±2- At larger s the effect is weak- 
ened because the total number of contributing subbands 
increases while the relative contribution of a particular 
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FIG. 2: (Color online) Spatial profile of the modulus of 
^n™(0, z) for subbands (0,0) (a) and (0, ±1) (b). The data 
for /x//kJ|| = 1.8, 2.2 and 2.8 are given in each panel and cal- 
culated at zero temperature. 



resonant subband diminishes. Results for a = — 180nm 
in Fig. [TJd) exhibit a similar behavior. 



Squeezing of fermionic-pair "wave function" 



The arguments given in Sec. IIII Al demonstrate that 
the redistribution of the fermionic condensate over the 
available subbands has a significant effect on ^q- An- 
other contributing factor is the large variation in the 
subband pair size ^|^"''"-' when the bottom of the cor- 
responding subband crosses ^, as seen in Figs. [Ijb) and 
(d). A further insight is obtained by considering how 
^,„„(r, r') decays with increasing z ~ z' . This decay 
of the subband pair wave-function ^„„j(r,r') is con- 
trolled by its longitudinal component ipnrniz, z') defined 
by Eq. (|3]). Figure [2] shows \ipnm{0, z)\ as a function of z 
for subbands (n, m) ~ (0, 0) [panel (a)] and (0, ±1) [panel 
(b)] at s = 1.8, 2.2 and 2.8 (for a ^ -140 nm). As 
seen, ■!/'o,o(0,z) is a slowly decaying oscillatory function 
for all values of s, which is typical for loosely-bound 
Cooper pairs in a bulk superconductor. Contrary to this, 
ipo,±i{0,z) exhibits a crossover from the strongly local- 
ized (at s = 1.8) to the BCS weakly localized regime (at 
s = 2.8) when the system passes through the resonance 
associated with s = 2. At s = 1.8 the lower edge of 
subbands {n,m) = (0, ±1) is situated slightly above /i, 
whereas at s = 2.2 and s = 2.8 it is slightly below and far 
below n, respectively. The edge of subband (0, 0) is far 
below jj, for all given values of s. Notice that at s = 1.8 
subband (0, 0) makes a contribution of about 95% to the 
total condensate while subbands (0, ±1) yield only 5%. 
At s = 2.2 the contribution of subband (0, 0) decreases 
down to 30%, as opposed to the contribution of subbands 



(0, ±1) that increases up to about 70%. 

The results of Fig. [5] arc understood as follows. When 
the lower edge of a single-particle subband is far below 
/i, the ratio of the pair-interaction energy J22| to the lon- 
gitudinal kinetic energy in this subband is small, as ex- 
pected in conventional weak-coupling BCS theory. How- 
ever, when the subband edge approaches /i, the longi- 
tudinal kinetic energy is reduced, and the ratio strongly 
increases, which means that the pair-interaction prevails 
over the longitudinal motion. This forces the fermionic 
pairs in this subband to squeeze in the longitudinal di- 
rection [see Fig. [2Ib)]. Thus, the effect follows from a 
redistribution of the kinetic energy between the longitu- 
dinal and transverse degrees of freedom in the subband 
whose lower edge crosses /i. This drop in the fermionic- 
pair size is in fact an atypical example of the BCS-BEC 
crossover that takes place in a single subband. However, 
unlike previously discussed systems, here the crossover 
is driven by the QS effects, which prompted the term 
atypical. Qualitative difference between the partial con- 
densates associated with different subbands explains the 
failure of the BCS estimate ^o ex Hvp/i^) in the analysis 
of the results in Figs.jTJb) and (d). 

Since the classical paper by Cooper [23| it is well-known 
that the configuration of the phase space available for 
the scattering of time-reversed fermions plays a crucial 
role for the formation of condensed fermionic pairs. In- 
deed, only a strong enough attractive interaction between 
fermions with opposite spin in 3D is able to produce a 
two-body bound state in the vacuum. However, when 
the scattering of fermions is influenced by the presence 
of a filled Fermi sea, i.e., the available phase space is re- 
stricted by exclusion of the single-particle states inside 
the Fermi sea, we arrive at the Cooper instability result- 
ing in the formation of weakly bound in-medium pairs of 
fermions for arbitrary strength of the attractive interac- 
tion. Restricting the phase space by removing the filled 
Fermi sea, one actually removes long range contributions 
in the Cooper-pair wave function, which, say, "encour- 
age" fermions to form in-medium bound states. Our 
results show that the additional reconfiguration of the 
phase space, such that the band of single-particle states 
splits up into a series of lower-dimensional subbands, can 
further modify the scenario of pairing through the atyp- 
ical BCS-BEC crossover. 



C. Coherent mixture of BCS and BEC-like 
condensates 

The BCS-BEC crossover in a resonant subband also 
reveals itself in the properties of the aggregate conden- 
sate of the system. However, the effect depends on the 
relative contribution of a resonant subband. It is maxi- 
mal at the corresponding resonance and diminishes away 
from the resonance. This contribution similarly drops 
when the total number of relevant subbands increases, 
and thus, as for other quantities, the lowest resonances 
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FIG. 3: (Color online) Contour plots of the order parameter 
A(p, z) obtained for a = —140 nm (at T = 0) below the reso- 
nance at ij,/hiuj_ — 1.9 (a), in the vicinity of the resonance at 
fi/h(jj± — 2.1 (b) and far beyond it at ^/hu}± — 2.8 (c). 
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FIG. 4: (Color online) Contour plots of the single-particle 
density np{p,z) (given in units of 10^^ cm^'*) calculated for 
the same parameters as in Fig. [3la)-(c), respectively. 



(i.e., for s = 2 and 3) are the most visible in ^o- Ef- 
fects of the subband BCS-BEC crossover on the aggregate 
condensate can be estimated by calculating the quantity 
7 = kp^o HI. For the BCS system 7 > 1, i.e., the 
Cooper-pair size greatly exceeds the mean distances be- 
tween particles, and this leads to a considerable overlap 
between the fcrmionic pairs. For 7 <C 1 such an overlap is 
absent, and the system becomes a EEC of tightly-bound 
point-like molecules. The intermediate crossover region 



is reached when I/tt < "/ < 27r [24 



For 



1.9 i.e.. 



just before the resonance at s = 2, we obtain 7 w 10 and 
so the system is in the BCS regime. At s = 2.1 wc find 
7 PS 4, which corresponds to the intermediate regime 
of the BCS-BEC crossover. The size-dependent drops 
in ^0 become larger for smaller \a\ (on the BCS side of 
the Feshbach resonance) , see Fig. [T] The reason for this 
is twofold. First, at smaller \a\ the energy window for 
contributing subbands also decreases, which means that 
the relative contribution of resonant subbands is larger. 
Second, the difference between ^^i^"™-' jn neighboring sub- 
bands increases, as seen in Figs. [TJc) and (d), and this 
increases the magnitude of variations in ^0 ■ 

Variations in ^0 are accompanied by substantial 
changes in the spatial profile of the order-parameter 
A(r) = — g\l'(r, r). Figure [3] shows the contour plots 
of A{p, z) for a = — 140nm and T = 0, calculated (a) 
slightly below the resonance, at s = 1.9; (b) close to 
the resonance, at s = 2.1, and (c) above the resonance, 
at s = 2.8. In the first case the pair condensate is al- 
most uniformly distributed over the trap with two peaks 
at the edges of the condensate cloud. These peaks are 
typical for a confined BCS condensate and can be ex- 



plained by the presence of the turning points in the tra- 
jectories of particles with energy close to p.. When the 
resonance develops [Fig. [3l[b)], the spatial distribution of 
the condensate acquires a pronounced bimodal charac- 
ter with an additional sizeable peak around z ~ due 
to the contribution of resonant subbands (0, ±1). The 
bimodal character clearly indicates that the system be- 
comes a coherent mixture of two qualitatively different 
condensates: the first is associated with subband (0, 0) 
and has properties typical for the BCS system; the sec- 
ond is due to resonant subbands (0, ±1) and displays a 
typical BEC-like behavior (the formation of bosonic-like 
states distributed at the center of the trap) . 

Note that contrary to the order parameter, the spatial 
distribution of atoms does not exhibit such noticeable 
changes at resonances, as seen from the comparison of 
Fig. [3] with Fig. |31 where contour plots of the position 
dependent single-particle density np{r) are given. Due 



to the axial symmetry, we have np{r) 
zero temperature 



np{p, z), and at 



^p{p,z) = XI (^ 



nmj 



An 



M 



E, 



iirnj 



l^nrn(p,^)Px?(^), (6) 



where the absolute value of the eigenfunction "dnmip, v) 
does not depend on ip. The sum in Eq. ([5]) is convergent 
and the ultraviolet regularization is not required. 



IV. DISCUSSION 

Here wc discuss details of approximations involved in 
our study, including the Anderson solution to the BdG 
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FIG. 5: The averaged energy gap (A) (in units of n) versus 
fi/huj±: squares are results for a — — 120nni calculated with 
the Bruun-Heiselberg cut-off |A^| < /i; circles are the data 
from Fig. [TJa) . 



equations, the ultraviolet regularization, and treating the 
interatomic collisions as a 3D process. 

The Anderson solution. Our calculations are based on 
the BdG equations that can be, in principle, solved with- 
out additional approximations. However, for the sake of 
transparent physical interpretations, here we employ An- 
derson's recipe for an approximate solution to the BdG 
equations, see Eq. ([T]). We recall that the Anderson ap- 
proximation in corp orates only the pairing of the time- 
reversed states |14l | , and the power of this approximation 
is based on the fact that the interaction matrix elements 
involving the time-reversed states are most pronounced 
as compared to other pairing combinations. Corrections 
to the Anderson solution were found in the range of a 
few percent for quasi-lD superconducting condensate in 
metallic nanowires [y] , which allows one to expect a sim- 
ilar accuracy of the Anderson approximation for a cigar- 
shaped ultracold superfluid Fermi gas. 

Notice that the energy spacing between the single- 
particle levels is not directly related to the validity of 
the Anderson ansatz. However, it can be an additional 
factor improving accuracy of the ansatz. Indeed, when 
the energy spacing between single-particle levels exceeds 
the characteristic pairing energy, the pairing of non-time- 
reversed states is suppressed not only due to smaller cou- 
pling but, in addition, due to a large interlevel spac- 
ing between the non-time-reversed states. In this case 
the Anderson approximation becomes practically exact. 
However, when the interlevel energy spacing is small as 
compared to the pairing energy, this does not necessar- 
ily mean a breakdown of the Anderson approximation. 
The pairing of the non-time-reversed states will not be 
pronounced anyway, thanks to the classical argumenta- 
tion by Anderson [1J|. For instance, Anderson's recipe 
yields exact solution for bulk superconductors where the 
pairing energy is much larger than the interlevel energy 
spacing. 

The ultraviolet regularization. Notice that our simpli- 
fied 3D-like regularization differs from the rigorous pro- 
cedure reported by Bruun and coauthors [15| for spatially 



nonuniform systems by the absence of an additional cor- 
rective term. However, following arguments by Bruun 
himself and Heiselberg in [23] [see the discussion just 
above Eq. (15) in this reference], one can expect that the 
role of this corrective term is not significant. Moreover, 
the authors of Ref . [25| introduced a simple cut-off ] A,^ ] < /i 
and argued that such a cut-off is "a first approximation to 
the more rigorous procedure of Bruun et al." for trapped 
Fermi gases. In Fig. [5] we show results for (A) calculated 
with this cut-off versus our data from Fig. [TJa) . Though 
(A) based on the cut-off procedure is slightly different as 
compared to the averaged gap calculated from Eq. ([2|), 
this difference practically disappears with a small shift in 
the scattering length, i.e., a -^ a + 20 nm, in the cut-off 
data. Thus, one can conclude that our simplified ultra- 
violet regularization, which becomes exact by increasing 
the number of subbands as the radial confinement is re- 
duced, does not affect our conclusions. 

Interatomic collisions. Several explanations about in- 
teratomic collisions in the cigar-shaped Fermi gas are also 
needed. For our trapping potential we have uju/lu± ~ 
0.01, and looking at this aspect ratio, one might get 
the impression that the two-particle scattering has an 
effectively- ID character here. However, this is not cor- 
rect. Though the character of the interatomic scattering 
is rather complex in quasi- ID systems, the choice of the 
3D pseudopotential is well justified for our range of pa- 
rameters. 

Based on the results for the binary atomic collisions in 
the quasi-2D tightly confined system |2a], one can con- 
sider that the scattering amplitude in the quasi- ID sys- 
tem is a function of the two important parameters ]a|/^_L 
and Esc/fiLUj_ (with Egc the scattering energy). Pro- 
nounced deviations from the ordinary 3D scattering can 
be expected [2^ when |a|/Z^ ^ 1 and Esc/hbj± ^ 1 (see 
alsoRefs.lllandlli). 

In our calculations l± varies from 0.35/im to 0.5/im 
when fi/huj± increases from 2 to 4. So, one finds that 
\a\/l± fa 0.3-0.4 for a = — 140nm. The relevant scat- 
tering energies in our problem can be roughly estimated 
as twice the chemical potential measured from the low- 
est single-particle energy, i.e.. Esc ^ 2(/i — haj±). Us- 
ing this estimate we find that for our parametric choice 
Esc/huj± ^ 2 (as // > 2huj±), which is directly related to 
the fact that our study is focused on the case of multiple 
contributing subbands. Thus, Esc/hujj_ is too large and 
a|/Zx is too small to favor pronounced ID-modifications 
to the interatomic collisions. We note that the choice 
of the 3D scattering length a = — 140nm and — 180nm 
is not crucial for our predictions of the atypical BCS- 
BEC crossover. We have performed additional calcula- 
tions and found even more pronounced variations of the 
fermionic-pair size at a = — lOOnm (see the discussion in 
the first paragraph of Sec. lIIIC"]) . 

It is also important to note that there is one excep- 
tion when the estimate of Egc in the previous paragraph 
does not hold. In a resonant subband, whose bottom 
is located in the vicinity of the chemical potential, the 



longitudinal motion of atoms is depleted. As a result, 
the relevant energies of the interband scattering in such 
a subband can be smaller than hLdj_. In this case mod- 
ifications to the pseudopotential could be pronounced, 
including the appearance of the confinement-induced Fes- 
hbach resonance [23, 123, if th^ ratio |a|//^ were large, 
i.e., \a\/l± > 1. However, for \a\/l± ^ 1 (this is our case) 
the ID-modifications are reduced to almost insignificant 
renormalization of a, see Ref. |271 

Thus, we can conclude that the ID-modifications to 
the 3D pseudopotential can be neglected for our choice 
of physical parameters. 



and BEC-like components, and the longitudinal pair size 
^0 associated with the aggregate condensate decreases. 
Radio-frequency spectroscopy can be used to detect the 
quantum-size driven squeezing of fermionic pairs |19|. 
Note that a similar many-body physics driven by the 
Quantum-Size oscillations can be expected for a pancake- 
shaped superfluid Fermi gas with only a few available 
transverse levels (the experimental realization of such a 
system, with observation of many subbands due to quan- 
tum confinement, was recently reported in [l2|). 
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